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O Abstract 

We consider correlation functions of the form (vacjOlvac)', where |vac) is the vacuum eigen- 
state of an infinite antiferromagnetic XXZ chain, |vac)' is the vacuum eigenstate of an infinite 
XXZ chain which is split in two, and O is a local operator. The Hamiltonian of the split chain 
(— I has no coupling between sites 1 and and has a staggered magnetic field at these two sites; it 

Q-i arises from a tensor product of left and right boundary transfer matrices. We find a simple, 

r~| exact expression for (vac | vac)' and an exact integral expression for general (vacjOlvac)' using 

the vertex operator approach. We compute the integral when O = af and find a conjectural 
^ expression that is analogous to the known formula for the XXZ spontaneous magnetisation 

and reduces to it when the magnetic field is zero. We show that correlation functions obey a 
boundary qKZ equation of a different level to the infinite XXZ chain with one boundary. 

> 

cn 

^ 1 Introduction 

The periodic XXZ chain is a quantum system with Hamiltonian 
^ N 

I: hxxz = -Iy: + -J+i^J + ^^l+i^l) ' (1-1) 

}_j where cr" denotes the Pauli matrix acting at the j'th site of the chain. We label our sites 



from right to left and identify the + I'th site with the 1st site of the chain. This simple model 
has been the subject of a huge amount of attention from the time of Bethe jlj until the present 
day (2). There are two main reasons for this interest. Firstly, the model is quantum integrable and 
has proved to be both a source and test-bed for a set of theoretical techniques applicable to such 
systems. These approaches include the Bethe Ansatz, quantum inverse scattering [s], the corner 
transfer matrix and Q operator Q , quantum field theory [5] , quantum groups [6|[7] and the vertex 
operator approach [8|[9j. Secondly, the XXZ model is of physical interest; it was developed as an 
idealised model of ferromagnetic systems and has recently found direct experimental realisation 



in a range of quasi-one-dimensional materials [10-12 . The Hamiltonian has also made a guest 



appearance in A/" = 4 super Yang Mills theory in the last decade 13 . 
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One of the key goals of all the theoretical work on the XXZ chain has been to exploit its 
quantum integrability in order to compute exact expressions for correlation functions. This goal 
was finally achieved for the infinite lattice in the 1990s. Two independent approaches were 
developed. The first was invented by the Kyoto group and deals with the infinite system with 
antiferromagnetic boundary conditions (8||9]. In this case, the Hamiltonian commutes with the 
action of the quantum affine algebra Uq (s[2 ) . This observation ultimately allowed the Kyoto group 
to produce an entirely representation theoretic expression for arbitrary correlation functions. 
We will refer to the techniques of the Kyoto group as the 'vertex operator approach'. The 
second approach, developed by the Lyon group, involved a direct Bethe-Ansatz computation of 
correlation functions for the finite-size system and a careful thermodynamic limit [l4 15 . 



In this paper, we study the XXZ model in a setting that we refer to as a fractured XXZ chain. 
We consider the eigenstates of two separate Hamiltonia. The first Hamiltonian is simply the 
infinite-size antiferromagnetic Hamiltonian 

that acts on the space with antiferromagnetic boundary conditions H 1 1 — or — I 1 h 

at plus and minus infinity (we use + and — to refer to two eigenstates of - see Section 2 
for a more precise definition of this space). This Hamiltonian is the one already dealt with by 
the Kyoto group and we simply recall the results that we need. In particular, labelling the two 
antiferromagnetic boundary conditions by z = and z = 1, we make use of the two corresponding 
lowest energy eigenstates |vac)(j) constructed in J9j. The second Hamiltonian we consider, which 
acts on the same space as above, is a split Hamiltonian 

H' = Hl + Hji, 

where the left and right Hamiltonia are 

1 

i<o 

With this Hamiltonian, there is split or fracture between the sites at position 1 and 0, with no 
nearest-neighbour interaction between them. There is also a staggered magnetic field h that acts 
at the boundaries of the fracture. Let us label the vacuum of H' by |vac)'^.^. This state will of 
course be simply a tensor product of the respective lowest eigenstates of Hl and Hr. 
The key results of this paper are exact expressions for 

(i)(vac|vac)(.), (1.2) 

and for all correlation functions of the form 

(i)(vac|C'|vac)(i), (1.3) 
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where O is any local operator. We obtain these expression by developing the vertex operator 
approach. 

There are several reasons why we look at the matrix elements (1.2) and (1.3) for this system. 
Firstly, they are 'natural' objects to consider in the language of the related 6-vertex model. 
We shall develop this point of view throughout the paper, but in brief, the matrix element 
(j)(vac|vac)j.^ is related to the partition function of the infinite lattice shown in Figurej2j As such, 
it can be expressed as a matrix element of two boundary, and two bulk corner transfer matrices 
representing the four quadrants of the lattice. Corner transfer matrices of both boundary and 
bulk types have been considered previously in the vertex operator approach [9, 16 , and it is very 
natural to put the two types together to form the fractured partition function of Figure [2| 

The second motivation for our interest in these fractured systems is their potential physical 
significance. One application is in the description of a local quantum quench in which a spin chain 
that consists of two separated domains for time t < is joined at time t = 0. Such local quenches 
have been investigated using conformal field theory techniques in 17-20 . In such a context, the 
squared modulus of the overlap (vac | vac)', which we compute exactly, is referred to as the fidelity 



and plays a role in understanding the dynamical behaviour of quantum entanglement 20 ,21 



A final reason for our interest in this fractured geometry is the potential connection with 
existing work on two-dimensional models on a wedge of angle a. Both massless and massive 



systems have been considered in such wedge geometries in 22 , 23 and [24| , and there exists the 
possibility of establishing a relation with our work when a — )• 27r. 

In this paper we construct the matrix element (j)(vac|C'|vac)|^j as a specialisation of a function 
P^^\(i, C2> • • • , Cn) defined as a matrix element of a product of vertex operators. We find that 
P^^\(i, C2, • • • , Cn) satisfies a boundary qKZ equation similar to the form found in 16 , but at 
a different level. We use a free-field realisation to derive a general multiple-integral expression 
for P(*)(Ci, C2) • ■ ■ ) Ca^)- Specialising to the case of the magnetisation, corresponding to the choice 
O = af, we integrate the resulting expression to produce a q-expansion to order q^^. Based on 
this expansion we make the following conjecture for the exact expression 

^2^n 



(0)( 



vac (Ti vac 



l + 2(l-r)^-^ 



where — 1 < r < 1 is related to the magnetic field /i > by 

h 



rq 



4n\ 



(1.4) 



1)1 +r 



Aq l-r 

Note that r = — 1 corresponds to the magnetic field /i = and r 
we have the special cases 



1 corresponds to /i = c«. Thus 



(0)( 



vac I CT]^ I vac 



'(0) 



n=l 



~ (-9^;?^)^' 



1, 



h 
h 



0, 



00 



The first expression coincides with the known bulk spontaneous magnetisation Q. We argue that 
this result and the h = 00 limit are as physically expected. We also make some observations 
about the special nature of the point r = 0. 
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The contents of the paper are as follows: in Section 2, we provide the definition of the model. 
In Section 3, we describe the realisation of the model within the vertex operator approach. We 
construct the vacuum and excited eigenstates for the two Hamiltonia H and H' and give a vertex 
operator expression for correlation functions. We then demonstrate how these same correlation 
function expressions arise from a corner transfer matrix approach. We go on to show that the 
correlation functions obey a version of the boundary qKZ equation. In Section 4, we make use of 
the free-field realisation to produce a general integral formula for correlation functions and show 



how its specialisation leads to the result (1.4). Finally, we summarise and make some concluding 
comments in Section 5. 



2 The Model 

In this section, we define our fractured XXZ/6- vertex model. All conventions closely follow those 



of the previous work on the vertex operator approach to the bulk [8|[9j and boundary 16 , 25 
XXZ models. Our starting point is to consider a vertex model defined in terms of the bulk and 
boundary Boltzmann weights R{C) ^-^id K{(; r) defined in [25] and reproduced in Appendix A. 
For convenience, we also define i^»(C) = K{(;r) and -R'o(C) = K{~Q~^C~^'if)- Let V be the 
two-dimensional space V = Cf + © Cw_ . The finite- lattice bulk and boundary transfer matrices 
are then operators — )■ defined by 



and 
where 



r/-(c) = Trvjno), 

no = i?oi(C)i?02(C)---^OJv(C) eEnd(yo©I^JV©•• 



(2.1) 



The boundary transfer matrix T^^{C,) [26] is represented by Figure 1, in which we make use of 
the pictorial conventions of Appendix A. We view this and other transfer matrices in this paper 
as acting in a downwards direction. 



1111 







c 



























Figure 1: The Transfer Matrix r|*"(C) 
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The two transfer matrices are related to bulk and boundary Hamiltonia if^*" and Hn^ by 



1 ^ 



(2.2) 



d 



Tjfin 



2q dC 
2 



logr^'"(C)lc=i + constant, 



^ N~l 



(2.3) 



,2 



where h 



Aq dC 
(g2 - 1) 1 + r 



"(01 (^=1 + constant, 



(2.4) 



4g 1 - r 

The partition function of interest to us is the infinite-lattice limit of the one shown in Figure 
2. We specify two different boundary conditions which we label by i = and i = 1. For i = 0, 
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Figure 2: The Partition Function of the Fractured Model 

the partition function is the sum over spin configurations of Figure [2] that are fixed at finite but 
arbitrarily large distances from the centre of the lattice to the ground state shown in Figure [3| 
For i = 1, the configurations are restricted by the corresponding ground-state with + -H- — . 

A key feature of the vertex operator approach is that it formulates such partition functions in 

(i) (i) 

terms of operators that act on the half-infinite left and right vector spaces n'£' and n'^' defined 

by 

= Span{---®t;,(2)^t;e(i)|e(n) = (-l)"+\n>0}, 
n'-^^ = Span{w^(o)®Ve(_i)0---k(«) = (-1)"^',«<0}. 

If we identify a dual vector w^- with f=p, then we can identify Ti^^ with 'H*j^^K The full horizontal 
space associated with the partition function of Figure 2 is 

= n^l^ ~ nf ® nf^ ~ Endin^l^). (2.5) 
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Figure 3: The i = Ground State 



There are then two different horizontal transfer matrices associated with the partition function 

associated 
associated 
>5 Tr{C) of 



represented by Figure 2: the bulk transfer matrix T{Q : J^^*-* — t- J^^-*^^*) of Figure 
with the lower half of Figure 2; and the transfer matrix T'(C) : J^^*) — t- J^*^*) of Figure 
with the upper half. We can clearly write T'{Q as the tensor product T'{Q = Tl{Q) 
the left and right transfer matrices Tl^r{C,) : T-L^j^^ '^^lr shown in Figure [ij 
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Figure 4: The Bulk Transfer Matrix r(C) 
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Figure 5: The Fracture Transfer Matrix r'(C) = Tl{0 ® Tr{Q 
Differentiating these transfer matrices, using the explicit form of R{C)i ^oiC) and K,{() given 



6 



in Appendix A, gives the associated infinite-lattice Hamiltonia H, H' : J^^*) — )■ J^^*) as follows: 

^ ■■= -^EKi-f + ^1+1^1 + ^+1-1) (2-6) 
^ \ogT{C)\(^=i + constant, (2.7) 



2q dC 

andH' := Hl + Hr, Hl,r : uf^^^ 

^^■■=-\T. (-^i-i + -J+i-l + ^-i+i-l) + (2 
1 



with . 

2 

i>i 



and ^il = "2 E (<-f-i + -i-l-i + ^-l-I-i) - ^-0, (2.9) 



where ff' = — -^—T' iC)\c=i + constant. (2.10) 
Aq dC ^ 

The first task in this paper is to diagonalise the Hamiltonia H and H'. Fortunately, this task 
has already been carried for both H and Hl by using the vertex operator approach to the bulk |9j 



and boundary 25 XXZ chains. It is then a simple matter to also diagonalise Hr by relating 
it to Hl', we describe the eigenstates in the next section. Denoting the unique, lowest-energy 
eigenstates of H and H' in the space J-"*-*-* by |vac)(j) and |vac)j^^, we can identify the value of the 
matrix element (j)(vac|vac)'^^.j with the partition function of Figure 2. 

The other objects of interest for us are correlation functions of the form 

(i)(vac|ii^j^ i.^, |vac)(.), 

£ (i) 

where E^, is the operator that acts at the j'th site of the lattice as 

El,{Va) = 5a,eVe'. (2.11) 

We identify this correlation function with the partition function of Figure [6] (shown for the case 
m = 2), in which we sum over only those configurations with the indicated spins fixed to the 
values {Emi • • • , ^2, El) and (e^, • • • , e'j, e'^). We again borrow the technology of the vertex operator 
approach to the bulk and boundary XXZ chains in order to compute these correlation functions. 



3 The Vertex Operator Approach 

In this section, we describe the realisation of the objects ^(0) ^'(C) and E^^j"^^ within the 
vertex operator approach. We rely heavily on the previous papers [9 
there for a fuller explanation. 

3.1 Identification of Spaces and Operators 

It is possible to define an action of the quantum affine algebra C/g(s[2) on the space by making 
use of the infinite coproduct. The essence of the vertex operator approach is that with respect to 



,IT6 |25l and refer the reader 
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Figure 6: The Correlation Function 



(i) 

this action it is possible to identify with the infinite-dimensional irreducible highest- weight 
C/g(s[2) module V{Ai) - where Ao,Ai are the fundamental weights. Similarly, the vector space 
T-L^^ ~ can be upgraded to a dual C/g(s[2) module (see Chapter 7 of [o] for the meaning of 
the * in this case) and is identified with V{Ai)*. Thus, there is a J7q(sl2) identification 

^« ~ V{A^) V{Ai)* ~ End(y(Ai). (3.1) 



The right-hand isomorphism in (3.1) is given by identifying a(^b* G V{A)(S)V{Ai)* with axb*{.) G 
End(y(Ai)). 

The next step is an identification of the half-infinite lattice transfer matrix which is shown in 
Figure "^a) and is a map from H^^^ — )• T-L^^ ^\ This is identified in terms of a uniquely defined 



(a) 
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1 


1 
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1 


I 


I 





























c 













(d) 



Figure 7: The Four Half-Infinite Transfer Matrices 

Uq^sh) map 

$(C) : V{Ai) ^ y(Ai_i) v^, 

where is the principal evaluation module associated with the two-dimensional Uq{sl2) module 
V = Cw+ ® Cv- [9j. Components $±(C) of this 'vertex operator' <I>(C) are defined by 
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and are vector-space maps V{Ai) — )■ l/(Ai_i). The associated transpose map V{Ai)* — y(Ai_j)* 
is denoted by <I>£(C)*- The left half-infinite transfer matrix of Figure [7][^ a) is then identified with 
the operator <l>e(C) '■ V{Ai) — )■ V{Ai^i), whereas the right half-infinite transfer matrix of Figure 
^h) is identified as ^-s{CY ■ V{Ai)* V{Ai^i)* (both identifications are up to a normalisation 
factor g2 defined by Equation (4.3) below). 

A dual vertex operator $*(C) : V^Ai) (8) V( — )• y(Ai_j), its components $j_(C), and transpose 
<l*j-(C)* are defined similarly in jOj. The two remaining half-infinite transfer matrices of Figure 
gc) and (d) are then identified as $*(C) : V{Ai) y(Ai_i) and <^>*_,{CY : ^(A^)* ^ V{Ai_i)* 
respectively. In summary, the four half-infinite transfer matrices of Figure [7] are identified with 
the four vertex operators 

(a) g'^MC), (b) ^?5$_,(C)^ 
(c) g'^m), (d) 5^f*_,(C)*. 

The half- infinite transfer matrices of Figure [7| and the reflection matrices K,{Q and i^o(C) are 
the building blocks of both T(Q and T'{(). Given the identifications of the previous paragraph, 
it is then straightforward to read off the corresponding realisation of r(C) and T'((") of Figures |4] 
and [5] in the vertex operator approach. We have 

r(C) = gJ^MO^'^-siCy ■■ V{Ai)0V{A^r ^V{A,^^)0V{Ai.ir, (3.2) 

e 

T'iO = Tl{C)'^Tr{C) : V{A,)CSViA,)* ^V{A,)(^V{Aiy (3.3) 
where TUC) = 5 E ^*'(^"')^-'(0'5.(C), (3.4) 

and Tnic) = gY.<^>iAcyKi'm-e{cy. (3.5) 

The right transfer matrix Tr(C) can be conveniently rewritten in terms of Tl{C,) by making use 
of the property 

K,{C) = K,[-q-\-^), (3.6) 
and the fact that ^^(C) = ^-e{~Qr^Q)- In this way, we obtain 

TR{0=TL{-q-\-^Y. 



Furthermore, Tl[C,) defined by (3.4) is simply the same as the boundary transfer matrix of [25] , 
there denoted by Tb{C)- It is shown in |25j that TB{—q~^C~^) = Tb{C)- Hence, in summary we 
have 

T'(C) = Tb{C)<^Tb{CY, where Tb{C) = 9j2^*e'((~^)^'t'iC)MC)- (3.7) 
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3.2 Eigenstates 

The eigenstates of the bulk transfer matrix T(() were first found in the paper [s], and are presented 
in the current notation in ^ . The form of the eigenstates is neatest if they are written in the space 
End(y(Aj)) ~ V^(A) (g) V{Ai)*. The state |vac)(j) is particularly simply and is just proportional 
to {—q)^, where D is a C/g(s[2) derivation that defines a grading of V{Ai) given by D\Ai) = |Aj) 
on the highest weight state |Ai) G V{Ai) and D fj^fj^ ■ ■ ■ fjJ\Ai) = m fj-^fj^ ■ ■ ■ fj,JAi) on a 
descendant state produced by the action of J7g(s[2) generators /o or /i. The inner product of 
two states in End{V{Ai)) is defined in joj by {f,g) = Try(^/^.j{f o g)^ and the normalised vacuum 
state is specified by 

|vac)(i) = ^(-g)^, with x = T^V{A,){{-qf'') = (3.8) 

for which r(C)|vac)(i) = |vac)(i_i), and (i)(vac|vac)(j) = -Try(A^) ((-g)^^) = 1. 

X 

We use the standard infinite-product notation 

oo oo 

(a;6i,62,--- ,bN)oc = J] 11 (1 " « ^2' " " • ^7)- 

711=0 71JV=0 

Also, note that we use the term eigenstate loosely here: |vac)(j) is only really an eigenstate of 
T^(C) because of the i ^ I — i shift above. It is however a genuine eigenstate of the XXZ 
Hamiltonian H given in Equation (2.6). 

Vacuum eigenstates \i)B £ ^(Aj) of the operator Tb(C) in the space V{Ai) were constructed 
in the paper |25]. They were found by solving the eigenstate condition 



rB(C)K)B = A«(C)| 



where the function A^'^\Q is specified in Section 2 of 25 . This equation is solved by using the 
free-field realisation of the module V{Ai). The dual vacuum state B{i\ was similarly found by 
solving the condition 



B 



{ilTsiO = A('\C) B{i\. 



Pull details, some of which are restated in Appendix B of the current paper, can be found in [25] . 

In our fractured XXZ model, it is the operator T'{() = Tb{C) 'S' Tb{CY that we wish to 
diagonalise. The form of this expression characterises T'(Q as an operator on the space V{Ai) (8) 
V{Ai)*; the associated action on an element / S End(y(Aj)) is simply 

r'(c)(/) = rB(c)o/ors(C). 

Hence, we can immediately write down the normalised vacuum vector |vac)j^^ of T'(C) as 

|vac)'.i) = \i)BB{i\ G End{V{Ai)), 
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or equivalently as 

|vac)^ = \i)B B{i\ e V{Ai) ViAi)*. 

We then have 

r'(C)|vac)'(,,) = -±^TB{C)\i)BB{i\TB{C) = {A^'HC)?W^c)[.y 

and (i)'(vac|vac)'(.x = Trv(A^(\i) BB{i\i) BB{i\) = l- 

(bW/b) ^ ' 

The overlap between the bulk and fractured vacua |vac)(j) and |vac)|^^ is given by the following 
expression: 



(i)(vac|vac)(.) = i ^ Trv/(A^) = ^ | B{i\{-q) 

X''bW)b ^ ^ X'^bW)b 



In the next section we use the free-field realisation to compute this overlap as a function of the 
magnetic- field parameter r. 

Excited states of the bulk Hamiltonian K are constructed in the vertex operator picture 
in [9]. The extra piece of technology required is new 'type IF vertex operators defined as suitably 
normalised [/^(sb) intertwiners: 

: y(A,) ^ y(Ai_,), ^{i) ■■ ^ n ® ^(Al-^), 

whose components ^'±(0 are defined by 

The key feature of these new vertex operators is that they quasi-commute with the 'type I' vertex 
operators $±(C)- Namely, 

^em%'{i) = T{<^liWAme[0. (3.9) 

*.'(0^£(c) = r(c/e)^e(c)^.'(e), (3.10) 



where t(^z) is the function 



z {qz-'^;q'^)ooiq^z'- 



'00 



It then follows from ( |3^ and ( [3^ that for / G End{V{Ai)), we have 

T{c)mof) = r{c/om)nc){f). 

Thus all states of the form 

€,(6) • • • €.™(6»^)|vac)(,) E End(y(A,)), 
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are eigenstates of T{() with eigenvalue t(C/^i) • • • r(C/^2m)- These states are conjectured to span 
the space ~ End(y(Ai)) in (o]. 



In the boundary case, it fohows from the definition of Tb(C) given by Equation (3.7) and from 



Equation (3.9) that 



TsiC) rsA^l) • • • re,J^2m)\i)B = liC, 6, " " " , C2n.W'HC) €,(6) " " " C™^™)!^)^, (3.11) 

where ^{C;£,i,-" ,S,2m) = Yl^i^iTiCCn^) '^iC~^S,n^)- III the context of the fractured model, in 



addition to (3.11), we also have 



^{i\A,ei^2i) ■ ■ ■ AA^i) TsiC) = 7(C; 6, • • • , ^2£)A«(C) B{i\A,M2e) ■ ■ ■ AA^i)- 



Thus a general excited state is of the form 
1 



B{i\i)B 

with T' {(,) eigenvalue given by 

7(C;6,--- ,6m)7(C;Cl,--- ,^2^) (A«(C))'. 



3.3 Correlation Functions 



The realisation of local operators E^.j'^^ (defined by Equation (2.11)) acting at the j'th site 



of the lattice in terms of type I vertex operators is given in Section 9.1 of [9]. The product 
j^e2(2)^£i(i) ijjgj^^i£g(j the following operator on V{Ki) (g) F(Aj)*: 



r U>_e;(-9"')«f-4(-g'') • • • ^-e'J-q-') ^.^(1) • • • ^e.(l)<I>.i(l) ® I 



(3.12) 



This realisation is the solution of the quantum inverse problem in the vertex operator picture. 
Let us define (for even) the following matrix element of vertex operators 

1 



^'^*^(Ci,C2,--- Xn) 



(i)(vac|vac)|.^ 



— (i)(vac|$(Ci)$(C2) • • • HCn) «> I|vac)'(^), 



(3.13) 



in which we are viewing the two vacuum states as a elements of V{^i) ® V{Ki)* and the product 
of vertex operators acts on the space V{Ki). Using the realisation of the vacua as elements of 
End(y(Aj)) given in the previous subsection yields the equivalent expression 



p(^)(Ci,C2,--- ,Cjv) 

or in components: 

^'^'^(Cl,C2,-- - ,C7v).i,.2,-,e„ 



1 



B(i\{-q)''HCimC2)---HQN)\i)B. 



(il(-g)^li) 



B 



B 



(3.14) 



B{i\{-q)''^eACl)^e,{C2)---^eACN)\i)B. 



B^\{-qr\^l 



(3.15) 
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Using the solution of the quantum inverse problem given by (3.12), we then have the following 
expression for correlation functions: 

1 



(i)(vac|vac)(.)^.^ 
with the choice Ci = C2 = ••■ 



■■ — Cm — —q ^, Cm+l — Cm+2 — •■■ — C2m — 1- (3.16) 

A general integral expression for P^^\(i, C2, • • • , Civ)ei,£2,-- ,£n given in Section [ijof this paper. 
3.4 The Corner Transfer Matrix Approach 



In Section 3.3, we have expressed correlation functions in terms of vertex operators via (3.16), 



where -P^*^ (Ci 1 C2 , • " " Xn) is defined by (3.13). Using the definition of |vac)(*) and |vac)'(*^ then 



gives expression (3.14) for P^^\(iX2, ■ ■ ■ Xn) in terms of the boundary states 1^)^ and the op- 



erator {—q)^. While the the argument leading to (3.16) is hopefully clear - at least when sup' 



plemented by the further details in [9,25 - the reader might still be lacking a simple, intuitive 
understanding of why correlation functions can be expressed via (3.15) and (3.16). Such an 
understanding is provided by considering corner transfer matrices. 

The idea of corner transfer matrices was introduced by Baxter, and their role in the vertex 



operator approach is described for bulk models in 91 and for boundary models in 16 . In brief. 



the 6-vertex model partition function which is associated with the Boltzmann sum over the 
infinite lattice with antiferromagnetic boundary conditions, and is represented by Figure [8j can 
be expressed as the trace of the product of four anticlockwise-acting corner transfer matrices 



(CTMs) A^^^^^{(), Agyy{C), Ag'^{C), A)!^'^{C) representing the four quadrants of the lattice thus: 



,(0 



7» 
^bulk 



n 



1111 



Figure 8: The Bulk Partition Function 



The CTMs act on the half-spaces H^l^ or Ti^^ as follows: 



/(») 



(i) 
L 



n 



L ' 



A%io : n 



L 



n 



(0 

R ' 



4^0 : n 



(i) 
R 



n 



(i) 
R ' 



4^0 : n 



(i) 
R 



n 



(i) 



13 



The remarkable property of corner transfer matrices discovered by Baxter is that in the infinite 



lattice limit we have ^(C) := ^sw(0 ~ C~^i where D : T-l^^ — t- Ti^ is an operator with non- 
negative integer eigenvalues and where ~ indicates equality up to a normalisation factor in this 
limit. We use the same notation D as for the t/q(s[2) derivation because we shall identify the 
two objects. The crossing symmetry of the 6- vertex R-matrix given by ( |A.2 ) then allows us 



to relate the other three CTMs to A{C) thus: A%{0 = SA{-q-^C~^), A^^>j^{C) = SA{C)S, 



^7vw(0 = '^C where S = ■ ■ - a^ a"^. Hence, we have 



"^bulk ~ 



The boundary partition function Z- 



boundary 



considered in 



16 



The modified CTMs associated with this model were considered in detail in 16 



25 



is represented by Figure 

The NW and 



1111 



} 



Figure 9: The Boundary Partition Function 

SW CTMs of Figure [9] were there labelled A^^y^^{(,l) and Agly{(,l), and the upper-boundary 
and dual lower-boundary states were denoted by \B; and (\ (the implicit ( dependence 

was suppressed in (I6]). In this boundary case we therefore have 

zi:Ldary = ^ (5; CI (C, i)4V (C 1) I i?; C) « • 

Identifying the boundary states as \i)B ~ ^ArV(C) 1)1^; C)^*"*) ~ C|^5vf(C; 1) then leads 

to 

^boundary ^ B{i\i)B- 

The partition function of the fractured model shown in Figure [2] is built by combining the 
CTMs and states of both the bulk and boundary models. It is apparent from Figure [2] that we 
have 

Zfracture = ^^^(5; Cl 4^ C, 1)4^(04^ (C)4W (C 1) I C)i^ 

where '•o^(-B; C| denotes the dual boundary state on the right of the fracture in Figure 2, A^^^{Q, 1) 
denotes the NE corner transfer matrix in the same figure, and \B;Q)^^ = \B-C)^\ Using the 
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crossing symmetry (A. 2) of the R-matrix and the relationship between Ko and K, given by (3.6), 



it then simple to show that we can identify 



C|4e(C, 1) = ^'HB; -q-^C'\AlU-q-'C\ l)S ~ B{i\S. 



Hence we have 



Z 



fracture 



(i3;Cl4;^(C,i)^«(-'z"k-')^«(C)4V(C,i)|i5;C)«~B(i|(-c?)^N 



B- 



Correlation functions associated with horizontal sites m, • • • , 2, 1 labelled as in Figure [6j can 
then be constructed in by inserting a product of 2m half-infinite transfer matrices of type (a) 
of Figure [7| beneath the NW CTM in any of the bulk, boundary or fractured scenarios above 
(the particular product of half-infinite transfer matrices required is always that of our expression 



(3.12) - a fact which again comes from the solution of the quantum inverse problem described 
in joj). Identifying these half- infinite transfer matrices with the vertex operators <^e(C)5 ^-nd H^^ 
with V{Ai), leads immediately to the general form for correlation functions (with even) 



G'«(Ci,C2,--- Xn) 



TVv(A,((-g)2^ 
1 



B 



B 

1 



Trv(A,((-'z)'''^(Ci)^(C2)---$(Civ)), 

B(i|$e,(Cl)'J>..(C2)---^e^(Civ)K>B, 



B{i\{-q)''\i)B 

found in the bulk, boundary and fractured models respectively. 



B(i|(-g)^$.,(Cl)^..(C2) • • • '^eACN)\i)B, 



(3.17) 
(3.18) 
(3.19) 



3.5 The Boundary qKZ Equation 

Boundary qKZ equations first appeared in the work of Cherednik |27| and were realised in a 



boundary lattice model in 16 . The two functions (3.18) and (3.19) obey boundary qKZ equations 



of different levels (we use the level terminology of |28[|29| ). This follows from the the following 
properties: 



K{cmc)\i)B 

PRiCl/C2)Hzi)^C2) 



A«(C;r)c/>(r^)|i)B, 
A(^){-q-\;r)B{imq-'C), 
A^'Hq-\;r)B{m-qfHq-\), 
HC2MC1), 



(3.20) 
(3.21) 
(3.22) 
(3.23) 



where the matrix KX) is defined by X) — ^-e'X)- The first two equalities (3.20) and (3.21) 
follow from the defining properties of \i)B and B{i\- 

TBX)\i)B = A^'HC)\i)B, and ^(^1 Tb(C) =B (i| A»(C; r) 

together with the following properties of vertex operators [9] : 



6e,e'l, $:(C) = ^-e{-q-\). 
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Equation (3.22) then follows from (3.21) supplemented by the vertex operator property {—q)^^{C) '■ 
^{~qC){~q)^ ■ The final relation (3.23) is the standard intertwining relation for vertex operators 
where P is the permutation operator P{vs-^ (8 ^^£2) — ^£2 

The boundary qKZ relations are simply the relations for the functions G(*)(Ci, C2; ■ " " ) Ca^) or 
P^'HCi, C2, • • • , Cn) above that follow from ( |3.20[ )-( |3.23[ ). In the case i = we have A(o)(C; r) = 1 
[25| and hence obtain 

G^^HCi, ■■■ , Cj-1) Q'~^Ci, Cj+1) • • • ; Cn) = 
R,,,.l{q-^Cj/Cj-l) ■ ■ ■ Rj,i{q-^Cj/Ci)Kj{-q~'Q 



X 
X 



^ij(CiCj) • • • Rj-i,j{Cj-iCj)Rj+i,j{Cj+iCj) • • • RnjiCnCj) 
K,{Cj)Rj,N{Cj/CN) ■ ■ ■ ii,,,+i(0/0+i)G'(°)(Ci, C2, ■■■ ,Cn), 



(3.24) 



and 



X 
X 



pW(Ci,-- - Xj-i,q~%Xj+i,--- Xn) = 
Rjd^iiQ'%/Cj-i) ■ ■ ■ RjA{q~%/Ci)Kj{q-\j) 

RljiClCj) ■ ■ ■ -^i-lj(Ci-lO)-Rj+lj(0+lCj) • • • RnjiCnCj) 

Kj{Cj)R,AQKN) ■ ■ ■ Rj,j+i{Cj/Cj+i)P^'HCiX2, • • • , Ctv). 



(3.25) 



Such generalised boundary qKZ equations were considered in [28]. In the language of [28] a 
boundary qKZ equation for a function ^'(Ci,C2j • • • Xn) is characterised as being of type (r, s) 
if we have ^(Ci,C2, • • • , Ca^) related to both ^(Ci, C2, • • • , r^s^/CN) and ^{A /(i, C2, • • • , Cn), as 
a consequence of which ^'(Ci, C2> • " " ) Ca^) is then related to ^'(Cij C2, • ' ' 1 '^'^Cj, Cn)^ Hence, we 
have (r, s) = {q~^, q^) for G^^\ and (r, s) = {q~^, q^) for P^^\ The 'level' £ of the boundary qKZ 
is then identified in I28j as s = (with k = 2 for the s/2 case) in analogy with the bulk 



qKZ equation of Frenkel and Reshetikhin |30|. Hence, in this language. Equation (3.24) (which 



appeared in [16]) is level 0, and (3.25) is level 2 



4 An Integral Formula for Correlation Functions 

In this section, we present a general integral formula for P^'^\(^i,C2, ■ ■ ■ ,Cn) which we then spe- 
cialise to give the fracture magnetisation 

c,-'\(vac|cjf |vac)/.x / ,^ ,^ \ 

I J^'^ = 5 p«(-5-S 1)-+ - P^H-<i-\ 1)+- • 

(i)(vac|vac)(.) V / 



We make use of the free-field realisation described in detail in the paper 25 and the book [9j. 
Most calculational details of the current work are relegated to our Appendix B to which to refer 
as necessary. 

^The slight difference in notation compared with f28l arises because z of 28 is equal to our 
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4.1 The Matrix Element (j)(vac|vac)'^^^ 

The boundary states are given in terms of highest-weight states |Aj) by \i)B = e^'|Aj) and 
B{i\ = (Aj|e'^% where Fi and Gi have the following quadratic form in terms of the q-bosonic 
oscillators a„ (n G Z/{0}): 



^ oo 

E 



t [2n][n] 



0, 



n=l 



[2n]N 



n=l 



The coefficients an, Pn ,^n)5n are functions of the magnetic-field parameter r entering the ma- 



trices K{C,,r) of Equation (A.l). They were defined in |25j and are reproduced in our Appendix 



B. 



B 



We show in Appendix B that we have 

(^2 10 fi fi\2 ^^2 6. „4 i 

\D\n.\ V 1 1 y 1 y Joo V y I y > y 7oo 



(0|(-g)^|0) 



B 



)2 

/oo 



B 



(0|0) 



(g2r2;g8)oo(g6;q8)^ 



and that and are given by the same expressions with the substitution 

r — )• r^^. Hence we obtain 

1 



(o)(vac|vac)(o) 



. , , , -B{^\{-q)''\i)B 



X^B{i\i) 



B 



{q^\q 



4\ 2 



°° (r2g4. ^8^ ^8)^(^2^12. ^8^ ^8)^ (^2^4. ^4^ ^8)^ (^10. ^8^ ^8)^ 



(4.1) 



The expression for (]^)(vac|vac)'^-|^^ is again given by the substitution r — )• r~^. 

The matrix element (o)(vac|vac)|Q^ represents the overlap between the bulk and fractured 
vacua, and it is of interest to see how this overlap varies as a function of the magnetic field h 
given by (2.4), for < /i < oo corresponding to — 1 < r < 1. The function (o)(vac|vac)|Qj is clearly 
a symmetric function of r with a maximum at r = 0. The value of the fracture magnetic field at 
this point is h = hmv ■= *''^4~^'* which corresponds to the Uq{sl2) invariant point 



31 



32 



of the 



finite boundary Hamiltonian (2.3). It is interesting that the vacua of H and H' are most similar 
at this very special value hinv of the magnetic field. The quantity |(o)(vac|vac)'^Q^p is called the 
fidelity in the language of quantum information theory, and a graph of the fidelity against h is 
shown in Figure 10 for two different values of the Hamiltonian anisotropy A = (g -|- q~^)/2. The 
fidelity increases with increasing |A| as can be seen. 



4.2 The General Integral Expression 

In Appendix B, we consider the free- field realisation of the general, N even, correlation function 

P^'\Cl,C2r-- ,CN)er,e„-,er, = rw \di A B (^| (-g)^1>e, (Cl)^^^ (C2) • • • (C^) 
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Figure 10: The Matrix Element |(o)(vac|vac)/Qj as a Function of Magnetic Field h 



(0)1 

for which we derive an integral expression. Letting A = < j < N,ej = +1}, and defining 



Zj = C|, we find 



N'^/A+iN/2- Yl a 



n 

j<k 



{q'^Zk/zj;q'^)c 



N l+e- 
l\{Wa- Wb){Wa- q'^Wb) 

^Tlh^l-' uT -2 ^nr (4.2) 



where l'^'^\{zj},{'Wa}) is given by Equation (B.4). The closed contour Ca is an anticlockwise 



one in the complex plane chosen such that the poles at q'^Zj lie inside and the poles at q'^Zj lie 
outside the contour. Further restrictions on the contour associated with the poles arising from 
l'^^\{zj}, {wa}) are given in Appendix B. Correlation functions are then given by the specialisa- 



tion corresponding to Equation (3.16) where the constant g is defined in j9j as 

4.3 The Fracture Magnetisation 

We consider the magnetisation 

m(vac|crf |vac)',.^ / ... 



(4.3) 



(j)(vac|vac)'. 



(i) 



Specialising formula (4.2) to = 2, with z = C, , we find 



dw 



gP'^H-Q-'CC)-^ = {q'zyzil - qr i ^ , (4.4) 

/(^C^) 2tt ^J —1 [w — z){w — q'^z)[w — q^z) 
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where /'» = F'(^ {q'z';q'U{qVz';q'U{q';q')M';Q')l&A'i'^') 

1 

^ 

(q'^zw; q^)oo{q^/{zw); q^)oo{q'^'^z/w; q^)oo{q'^w/z; q^)oo 

1 



{q'^zw; q^)oo{q^/{zw)] q^)oo{q^^z/w; q^)oo{q^w/z; q^)^ ' 

{q'^rz; q^)oo{q^r / z; q^)oo {q^rw; g^)oo(g^r/-w; g^)oo 
{q^rz;q^)oo{q^r/z]q^)oo {rw;q^)oc{q'^r/w;q^)oc 
(l/(rz); q^)ooiq^z/r; g^)oo {q^w/r; q^Uiq^jrw); g^)oo 
(rz); g8)^(g4^/^. ^8)^ [q'^w/r; q^)oo{q'^/{rw);q^)oo ' 



and 



The contour Cli^ is fixed fohowing the prescription discussed in the previous subsection and in 
further detail in Appendix B. Let us use the notation yi < w < y2 to indicate a w contour which 
passes outside the yi pole and inside the 2/2 pole. Then, to be explicit, the contour is given by 
the following requirements for all n > 0: 

Ci°| : q'^^'^/z, g6+«7z, gi2+8n^^ ^4^ ^ 

w < q-^yr, q-^-^'yz, q'^^^'^/z, q^-^'^z, q'^'^'^z, z, zq^ , 

CW . g2+8n/^^ ^4+8n/^^ q^+^^/z, q'^'+'^z, q'^+'^^z, q'z < w, 

w < q-^-^'^r, g-^-^'^/z, g-^-^"/z, g*-^"z, g-^'^^z, z, zg^. (4.5) 

Similarly, we find 

9P^H-q-\.0^- = -{q'zrz{l-q'f i '''~\. ^/'« 

JqM_ 27rv-l {w - z){w - q^z){w - q^z) 

(i) (i) 9 

where C_j__ is the same as C__^ except that the pole at q z now lies inside the contour. 

A useful check of our formalism is provided by computing the matrix element of the identity 
operator I^^) = E^^^^ + EZ^^^ . We have 



(j) (vac|l(^)|vac>(i) 
(i)(vac|vac)|.-, 



9 



(p»(-g-i,l)_+ + P«(-g-\l)+_). 



We see from the above discussion that 

g {p(^\-q-\l)^+ + p('\-q-\l)+\ = -{q^zYzil - q^fRes^^,,2 ' 



w 



{w — z){w — q^z){w — q^z) 



We find that this expression has the required value of 1 for both i = and i = 1. 

It remains to compute P^^\—q~^,l)^ . The obvious way to carry out this integral is by 

(i) 

a brutal summation of the infinite sets of poles lying inside C_j__. Another way is to recall the 
origin of the w integrals in the expression for ^'^*^(Ci) C2) • " " ,CN)ei,e2,--,epf- they arise because 
the vertex operator <I>+(C) is expressed as the (-l)th coefficient of a Laurent expansion in w (see 
Chapter 8 of (9|). At any given order in q, there are only a finite number of contributions to 
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this coefficient that come from expanding the numerator and denominator of the integrand of 
ei,e2,--- ,ejv (t&king care to respect the analyticity requirements expressed in the 
specified contour). Hence, the q-expansion of the correlation function is relatively easy to compute 

from the integral expression. In this way, we have computed up to order q^^. 

The result is consistent with the following conjectural form: 

,p(0)(-,-M;r),„ = (l-.)^J^, 

n=0 y / 

and also with the required symmetry 

P»(-g-\l;r)+_ =p(i-)(-<7-i,l;r"i)_+ (4.6) 

(we are now showing the implicit r dependence of the correlation function). The sum can be 
represented as a basic hypergeometric function (see, for example. Chapter 17 of |33| ) in the 
following way: 

Hence we obtain the magnetisation 

-MW(r) = g(^-P^'\-q-^C,C;r)^++P^'\-q~'C,C;r)+-)=-l + 2gP^^\-q-\l;r)+_ 

which can also be re-expressed at 



l + 2(l-r)E71^ (4.7) 



It also follows from the symmetry (4.6) that we have 



The form ( |4.7[ ) is similar, but different, to the magnetisation at site 1 of the pure boundary 
model of Figure [9] considered in [25] : 



<i» = 1+2(1 (4.^ 

n=l ^ y / 



The three special cases r = —1, 0, 1 corresponding to /i = 0, /ijn^, oo are of separate interest: 
When r = —1 {h = 0), we can sum the expression (|4.7[) to obtain 



-M(0)(r = -1) = 7^44f^- 



loo 



This coincides with the expression for the spontaneous magnetisation of the bulk model Q (the 
minus sign is just due to the choice of convention for the i = ground state of Figure Isl). It seems 
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physically reasonable that the magnetisation at position 1 in Figure |6] is simply the same as that 
of the zero-field bulk model when the fracture magnetic field is zero. 

• When r = {h = hinv), we have 



-M(°)(r = 0) 



which coincides with the magnetisation (4.8) of the boundary model at this value of the magnetic 



field. At this point the overlap (o)(vac|vac)'^Q^ is a maximum and thus the bulk system is closest to 
the tensor product of right and left boundary systems. It is therefore reasonable that the fracture 
magnetisation at site 1 should coincide with magnetisation at the corresponding site of the left 
boundary system shown in Figure [9j 

• When r = 1 (/i = oo), we have M(°)(r = 1) = -1. A gain, this is as expected, since the presence 
of the term haf in H' will force the magnetisation at site 1 to — 1 as /i — )• oo. This is also the 
value of the boundary magnetisation M^^^^{r) as /i — )• oo. 

In order to summarise the behaviour, we fix A = —2 and plot the two magnetisations M^^\r) 
(fracture mag) and M'^2i^^{r) (bound mag) as a function of h{r) in Figure [ll| We also show 
minus the value of the zero-field bulk spontaneous magnetisation (spon mag) for comparison. 
The coincidence of M'^^\r) with the spontaneous magnetisation at /i = 0, and with M^2i,^^{r) at 
hiriv and /i — )• oo is evident. 




spon mag 

bound mag 

- fracture mas 



Figure 11: The Spontaneous Magnetisation (with added - sign), and the Fracture Magnetisation 
M^^\r) and Boundary Magnetisation -^bo«nd(^) ^^^^ 1 as a Function of Magnetic Field h (at 
A = -2) 
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5 Summary and Discussion 



In this paper, we have addressed the novel problem of computing XXZ correlation functions in the 
fractured geometry of Figure [6j We have expressed these correlation functions in terms of matrix 



elements P^^\(i, (2, - ' ' > Cn) defined by (3.14) that obey a level 2 boundary qKZ equation. Making 



use of the free-field realisation, we have derived the general multiple-integral formula (4.2). We 
have specialised this formula to the case relevant to the magnetisation at site 1, and by carrying 
out a high-order q expansion of the integral have arrived at the conjectural form given by Equation 



(4.7). This form is similar, but different, to the magnetisation (4.8) at the corresponding position 
in the boundary geometry of Figure [oj The latter magnetisation was computed in [25| , in which 
case the corresponding integral was carried out exactly by summing residues. However, while the 
final magnetisation expressions are similar, the integrand in our fractured case is considerably 



more complicated that in the pure boundary case of 25 . Ultimately, this fact is linked to the 



additional {—q)^ factor that distinguishes the boundary and fracture correlation functions given 



by (3.18) and (3.19). The origin of this extra factor lies in the presence of two extra two corner 
transfer matrices A^g^^{()Ag^^{Q in the fractured partition function of Figure [2] when compared 
to the boundary partition function of Figure |9| This additional factor {—q)^ factor is also 
responsible for the different level of the boundary qKZ equations in the two cases. 

The corner magnetisation in a wedge geometry has been considered before. Conformal field 
theories in such a geometry were analysed by a conformal mapping from the upper-half plane [23] . 
For the Ising model on a triangular lattice, a finite-size scaling approach was taken in |24| , while 
the Ising model on a cone was considered in the same paper by taking a trace over a variable 
number of corner transfer matrices. It would be interesting to try to connect our magnetisation 
results with these existing calculations in the limit as the wedge angle approaches 2it. 

Almost as a bi-product of our analysis we have calculated the overlap (j)(vac|vac)j.^. As we 
have mentioned in the introduction, our fractured geometry, and in particular this overlap, have 
recently become objects of interest in the literature on local quantum quenches and their use as 



a probe of the dynamics of quantum entanglement [18 20 21 . So far, most results in this field 
have been obtained using conformal field theory techniques, and we hope and anticipate that our 
exact massive lattice model results will be valuable in this arena. 

Finally, the vertex operator approach to boundary problems has been generalised in a number 
of directions (34] - |37| , and there will be corresponding generalisations of the analysis of fractured 
models that we have developed here. In particular, the analysis should be generalisable to ABF 
models p4^ and higher spin models. The former generalisation might be particularly useful in 
establishing connections with conformal field theory results for wedge geometries. 
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A Bulk and Boundary Weights 

We make use of the following bulk and boundary weights |25|: 



(l 



V 



K{C;r) 



'l-rC^ 







f{C;r) \ 1, 



(A.l) 



where 



iq^C~^;g*)oo{gK^;g^ 

Components are defined by 



The matrix R{0 obeys the Yang-Baxter, crossing and unitarity relations, and i^(C; obeys the 
boundary Yang-Baxter, boundary unitarity, and boundary crossing relations shown in ^25j. In 
this paper, the property that is used most is the crossing symmetry of the R-matrix: 



RlKl{C) = R 



1 



(A.2) 



We make the further definitions 

K,{0 = K{C;r), and i^o(C) = K{-g~\-^;r). (A.3) 
The graphical representation of these matrices used throughout the current paper is that 



shown in Figure 12 



^43(Ci/C2) 



4 







Cl 








e'l 



£2, K,l'iO- 



Figure 12: Graphical Representation of the Boltzmann Weights 
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B The Integral Formula for Correlation Functions 

In this appendix, we give the details of the free-field calculation of b and of 

^'(*^(Cl,C2,--- ,CN)e„e„-,e^ = B (^| (-g)^"^.! (Cl)^e. (C2) • • • (Civ) I^B, (B.l) 

B{i\[-q) \i}b 

when N is even. We closely follow the approach and notational conventions of |25[ in which 
B{i\i)B and ,.|., siil^ei (Ci)*^'£2 (C2) • • • ^e]v(Civ)K) B are calculated; the key difference is that we 



replace B{i\ of 25 with g)-^. 



The states \i)B and are expressed in the free-field realisation of [25| as 
\i)B = e^'\Ai), B(i| = (Aile'^S 

1 ^ nan 2 , \ ^ n(i) ^ 1 "Tn 2 , \ ^ 

where = ^ L Mn}"'" + L/^i^"- = 2 Mn]'^'^ + ^ 

n=l J ra=l n=l ^ ^ n=l 

Here, the coefficients are given by 

an = V", 7n = 

(0) _ ^ g5n/2(l_gn) g7n/2^n ^ ^ ^-3n/2 _ ^n) ^-5n/2^; 



2n] [2n] ' " " [2n] [2n\ ' 

g3n/2^-n. (1) _ g-3n/2(;L _ ^n) g-n/2^-r, 

2^^^ ^ [2n] ' " ~ " ^^^^ ^ [2n] 



where 



g« _ I 1 if n is even, 

[«] ■= ZT' 8-^^ •= { 

Q-<1 if n is odd. 



Given the free-field realisation properties {—q) ^an{—q)^ = q^^dn and (Aj|(— g)^ = (Ajj, it 
follows that we have 



B{i\{-qr = (Aden's where G', = ^^^a2+j;5:«a„, 7; = g^"^^, = g^n^W. 

^-^ \Zn\\n\ ^-^ 

n=l ^ -I n=l 



Now we closely follow the procedure of Section 4 and Appendix C of i25!j with the change of 
In 7n and b'^n — >• bn^ . In this way, we obtain 



B 



{^{-qn)B = f[ '—r exp £ -J^^Ainm' + Pii' + ^'c^n) 



(B.2) 

The sums here can be written as infinite products using the identity 



00 



exp 



Carrying this procedure out and making use of standard q-product identities, as well as the 
less-standard, but easily-derived, identity 



(o; 6, 6)00 = (a6; 6, }?)\^ (a; 6^)00, 
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we find 



B{U\[-q) \U)b - 



and that is given by the same formula with the substitution r ^ r ^. Using the 

expressions 

(A'i^X ,,,,, {q^r~^;q'')oo 



of 



25 



, and the definition (3.8) of Xi we obtain the result for (j)(vac|vac)'^^^ given by Equation 



(4.1). 



Let us carry on to calculate the correlation function (B.l). Defining the set A by 

^ = {j|l<i< A^,ej = +1}, 
we again follow the method of Section 4 of |25] to obtairj^ 



i=l j<k 



\[{Wa-Wb){Wa- q^Wb) 



X 



27rV^ n {^3 - q-^Wa) n (^a " q^Zj) 

j<a a<j 



where the contour Ca is given by q^zj < Wa < (f'Zj (recall - the notation A < Wa < B means that 
the pole at A lies inside, and pole at B lies outside the anticlockwise Wa contour) and 

+ \anY^ + + 7;/3i*))^n - (/3« + an5;«)y„|) , (B.3) 



with 



7n/2 n/2 -5n/2 ^ n/2 



[2n]'^^j" [n] [2n] [n] 

Computing the sums as before we find 

l'(^\{z,},{Wa})=F'(^H{z,},{Wa}) 

yr iq''z,Zk;q\q')oo (gV(^,^fc); g^ g^)oo yr iq''> z, / z^; q\ q^'U yr jq'^ z] ; q^U iqVz];q\q^U 
""/i iq''^jZk;q\q^)oo {qyiz,Zky,q\q^)oo (q'^z, / z^; q\ q^)oo \^ (q^^ z] ; q^ q^U (q'/zj^q^q^U 

n (g6^,u;.;g8)oo(gV(^.^«a); 98)00^2^,/^,; 58)^(g6^^/^^..g8)^ n(^'^- ('^'^^i 9')oo (B.4) 

J, a a 

X JJ(g^^^aWfe; q^)oo{q^WaWb] q^)oo{q^/{WaWb); q^)oo{q^/ iWaWb);q^)oo Yliq^'^a/Wb] q^)oo{q^^'Wa/Wb; q^)oo, 
a<b a,b 



^There is a minor typographical error in (4.4) of 25 in which an additional factor is present on the left-hand 
side. However, this error is not carried through and the final expression (4.8) in ,25, is correct. 



25 



where 



F''^'H{z,},{Wa}) = 

J-.J- {q^rzj;q'^,q'^)oo {q^r / Zj;q'^,q'^)oo^} {rWa, q^)oo{q^r /wa^q^) oo 



f'W({z,},{^4) 

TT {'^1 {rzj)\q* ,q^)oo{q^ / {rzj)-q^ ,q^)^ (g^^j/r; g^ g^)^ -[-[- (q'^Wg/r; q^)oo{q^ / (rWq); q^ 

(gVK); g^ 'Z')oo (q'zj/r; q\ q^U V i<l^wJr; q'UiqVirWa); q^)oo 



oo 



The full contour Ca is now fixed by the previous requirement q'^Zj < Wg < q^zj supplemented 
by the requirement that the contour be consistent with the region of validity of each of the sums 
contributing to the last expression for /'(^^({zj}, {wg}). Hence, if a term (wg/A; q^)oo appears in 
the denominator of /'(^^({zj}, {wg}), we insist that Wg < Aq~^^ for n > 0; if a term {B/wa;q^)cx> 
appears we insist that Bq^^ < Wg for n > 0. 
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